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Abstract 

We study the geometric and electronic structures of silicene monolayer using density functional 
theory based calculations. The electronic structures of silicene show that it is a semi-metal and the 
charge carriers in silicene behave like massless Dirac-Fermions since it possesses linear dispersion 
around Dirac point. Our results show that the band gap in silicene monolayer can be opened up 
at Fermi level due to an external electric field by breaking the inversion symmetry. The presence 
of buckling in geometric structure of silicene plays an important role in breaking the inversion 
symmetry. We also show that the band gap varies linearly with the strength of external electric 
field. Further, the value of band gap can be tuned over a wide range. 
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I. INTRODUCTION 



There has been a lot of interest in sihcene since it shows properties similar to those of 
graphenei"-. For example, the theoretical studies on silicene show that the charge carriers in 
this two-dimensional material behave like massless Dirac-Fermions due to presence of linear 
dispersion around Fermi energy at a symmetry point K in the reciprocal lattice^"^. Similar 
to its carbon counterpart - graphene, silicene is a potential candidate for applications in 
nanotechnology. The silicon based nanostructures have an important advantage that they 
are compatible with the existing semiconductor technology. Therefore, silicene and silicon 
nanoribbon have received much attention from both experimentalist and theoreticians^"—. 
Recently, silicene has been epitaxially grown on a close-packed silver surface Ag(lll)^. 
Though graphene possesses many novel properties, its applications in nanoelectronic de- 
vices are limited due to its zero band gap and hence it is difficult to control the electrical 
conductivity of graphene. However, band gap in graphene may be introduced by chemical 
doping but chemical doping is uncontrollable and incompatible with device applications. 
Hence, it is desirable to have band gap in materials in addition to their novel properties. 

In this letter, we demonstrate that silicene is one such a material in which the band gap 
can be opened up as well as varied over a wide range by applying external static electric 
field in a direction perpendicular to the plane of sheet. We use SIESTA package^^-iS for 
performing a fully self-consistent density functional theory (DFT) calculation by solving the 
standard Kohn-Sham (KS) equations. The KS orbitals are expanded using a linear combi- 
nation of pseudoatomic orbitals proposed by Sankey and Niklewski^S. All our calculations 
have been carried out by using triple- zeta basis set with polarization function. The standard 
norm conserving TrouUier-Martins pseudopotentials^ are utilized. For exchange-correlation 
potential generalized gradient approximation given by Perdew-Burke-Ernzerhof^^. has been 
used. A cutoff of 400 Ry is used for the grid integration to represent the charge density 
and the mesh of k-points for Brillouin zone integrations is chosen to be 45x45x1. The 
convergence criteria for energy in SCF cycles is chosen to be 10~^ eV. The structures are 
optimized by minimizing the forces on individual atoms (below 10~^ eV/ A). Super-cell 
geometry with a vacuum of 14 A in the direction perpendicular to the sheet of silicene is 
used so that the interaction between adjacent layers is negligible. 

The optimized geometry of silicene is shown in Fig. 1. The unitcell has two Si atoms (A 
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FIG. 1: The optimized structure of silicene monolayer, (a) Top and (b) side view of 5x5 supercell. 
The hnes in (a) represents the unitcell. The vertical distance between two Si atoms at sites A and 
B is represented by 'd'. 

and B) and the space group is P3ml. We observed that the minimum energy structure of 
the sihcene is low-buckled with the lattice constant of 3.903 A. The bond length and bond 
angles between the silicon atoms are 2.309 A and 115.4° respectively. The value of bond 
angle in silicene lies in between those of sp^ (120°) and sp^ (109.47°) hybridized structures. 
This clearly shows that the hybridization in silicene is not purely sp^ but a mixture of 
sp^and sp^. The buckling in silicene is due to weak tt - n bond that exists between Si atoms 
since the Si-Si distance is much larger as compared to that in graphene (C-C =1.42 A). The 
system increases its binding energy due to buckling by increasing the overlap between n and 
cr orbitals. The presence of buckling may also be explained by Jahn- Teller distortion. Our 
results on geometry match well with the previous calculations^'^iii"— . Here, we quantify the 
amount of buckling in terms of buckling length d which is defined as the vertical distance 
between atoms at sites A and B in the unitcell. The value of d in silicene is 0.5 A. In case 
of planar structure, as that of graphene, the value of d is zero. The importance of d in 
determining the electronic structure of silicene under the influence of external electric field 
will be described later. 



The band structure along high symmetry points in Brillouin zone and density of states 
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FIG. 2: (a) Band structure and (b) density of states for optimized structures of silicene. The 
energy of bands are with respect to Fermi level. 

(DOS) of silicene are shown in Fig .2. The figure clearly indicates the semi-metallic behavior 
of silicene since the value of DOS at Ei? is zero and the conduction and valence band touch 
each other only at the symmetry K point. The energy levels and the contribution of density 
of states just below and above Fermi levels are mainly due to tt and n* orbitals. The 
dispersion around K point near the Fermi level is linear (see insert in Fig. 2 (a)). The point 
in E-k diagram where the conduction and valence band touch each other at E^ is called 
the Dirac point. The presence of linear dispersion indicates that the charge carriers near 
Dirac point behave like massless Dirac Fermions since the dynamics of these carriers obey 
relativistic Dirac-like equation. The relativistic Dirac-like Hamiltonian which describes the 
electronic structure of silicene around the Dirac point, similar to that of graphene^^i^, can 
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be approximated as 
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where k and vp are momentum and Fermi velocity of charge carriers near Dirac point. The 
quantity A is the onsite energy difference between the Si atoms at sites A and B. Due to 
the presence of inversion symmetry, the onsite energy difference A becomes zero, which 
leads to the linear dispersion around the Dirac point linear i.e. E = ^hvpk. Until now, 
the treatment for description of electronic structures of silicene is exactly similar to that of 
graphene. The band gap in graphene-like structure can be opened up if one can break the 
inversion symmetry and hence the value of A becomes finite. Then, the dispersion around 
the Dirac point becomes 



In this case, the value of band gap opened is twice that of the onsite energy difference. 
It is possible to make A non-zero in graphene by chemical doping which cause the local 
environment around the sites A and B to be different. As mentioned earlier, tuning the 
value of A by chemical doping is a difficult task. In the above mentioned discussion, the 
effect of spin-orbit coupling (SOC) is not included and very small gap of the order of /ieV 
in graphene and about 1.55 meV in silicene may be opened up due to SOO^"^. However, 
presence of small gaps may be useful in studying fundamental properties like quantum spin 
Hall effect and will not be useful in operating or controlling nanoelectronic devices made 
up of graphene/silicene at ambient conditions. In order to utilize many novel properties of 
graphene-like materials in nanodevices, it is desirable that the materials shall possess value 
of band gap more than the thermal energy kT at room temperature (~ 25 meV ). 

In Fig. 3, we plot the spatial distribution of highest occupied and lowest unoccupied 
states with and without electric field. Both of these states lie at the symmetry K point 
in Brillouin zone. In the absence of external electric field, both the highest occupied and 
lowest unoccupied states ( Fig. 3(a) and 3(b)) have same energy and also possess inversion 
symmetry. Hence the spatial distribution of these states above and below the silicene sheet 
looks similar. However, application of external electric field breaks the inversion symmetry 
in these states. This makes the spatial distribution of charges above and below the silicence 
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sheet to be different which can be seen clearly from Fig. 3(c) and 3(d) respectively for 
highest occupied and lowest unoccupied states. This leads to the opening up of band gap 
in silicene. 
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From our calculations on electronic structure of silicene, we observed that a band gap 
can be opened up at the Fermi level due to the influence of external electric fleld. The 
electric fleld is applied in a direction perpendicular to the plane of silicene sheet. The 
calculated band structures of silicene around Dirac point for various strengths of external 
electric fleld are shown in Fig. 4. For example, the value of band gap opened with electric 
flled of 5 volts per nanometer (V nm~^) is 50.9 milli electron volts (meV). The band gap 
in silicene is opened due to breaking of inversion symmetry by the electric fleld since 
the potential seen by the atoms at the sites A and B are different. In this situation, the 
finite value of onsite energy difference arises due to the potential difference and hence, we 
can write A = a{VA — Vb) where a and Va (Vb) are proportionality constant and the 
potential seen by the atom at the site A(B) respectively. When the fleld applied between 
the layers is constant around the sheet, as in the present case, the potential difference 
becomes (Va — Vb) = Fd, where F is strength of electric fleld and d the distance between 
the atoms at the sites A and B along the fleld direction. The flnite value of d in silicene 
monolayer plays an important role in breaking symmetry in presence of electric fleld and 
cause opening of band gap. This result leads to an important advantage of silicene over 
graphene since there is no buckling in latter and hence it is not possible to open up a gap 
by applying exterr 




FIG. 5: The variation of band gap (E^), valence (E^,) and conduction (Ec) band edges of silicene 
with the strength of external external field applied perpendicular to the plane of silicene sheet. 



We also observed from our DFT calculations that band gap can be tuned over wide range 
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by varying strength of external electric field. The results of variation in band gap of silicene 
under influence of different strengths of electric field are summarized in given Fig. 5. It 
is seen from Fig. 5 that the band gap varies linearly with the external electric field. We 
also observed that both valence and conduction band edges move symmetrically away from 
Fermi level when the electric field is applied. The presence of linear relationship between 
the band gap and electric field can be explained as follow. As mentioned earlier, the value 
of band gap opened can be written as Eg = 2A = 2a{VA — Vb) = 2aFd= a'F, where a' 
is proportionality constant and it should be characteristic of the material. From Fig. 5, 
the value of the a' is found to be 10.14 (meV per V/nm). It is remarkable that in silicene 
monolayer, the band gap can be tuned by simply varying the strength of external electric 
field. Similar trends have also been observed for germanene monolayer and the results will 
be published elsewhere^^. 

In summary, we have carried out abinitio DFT calculations to study geometric and elec- 
tronic structure of silicene. Our results on the electronic structures of silicene monolayer 
show that it is a semi-metal and possesses linear dispersion around Dirac point similar to 
graphene and hence the charge carriers in silicene behave like massless Dirac-Fermions. We 
have shown that band gap in silicene monolayer can be opened up by applying external elec- 
tric field. The presence of buckling in optimized structure of silicene plays an important role 
in breaking inversion symmetry. We also observe that the band gap varies linearly with the 
strength of electric field. Furthermore, our results predict that the band gap produced can 
be more than the thermal energy by applying electric field strength of few V/nm and hence, 
there is a possibility of using silicene monolayer in nanodevice even at room temperature. 
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